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ABSTRACT 



Context. Temperature, surface gravity, and metallicitity are basic stellar atmospheric parameters necessary to characterize a star. There 
are several methods to derive these parameters and a comparison of their results often shows considerable discrepancies, even in the 
restricted group of solar-type FGK dwarfs. 

Aims. We want to check the differences in temperature between the standard spectroscopic technique based on iron lines and the 
Infrared Flux Method (IRFM). We aim to improve the description of the spectroscopic temperatures especially for the cooler stars 
where the differences between the two methods are higher, as presented in previous work. 

Methods. Our spectroscopic analysis is based on the iron excitation and ionization balance, assuming Kurucz model atmospheres in 
LTE. The abundance analysis is determined using the code MOOG. We optimize the line list using a cool star (HD 21749) with high 
resolution and high signal-to-noise spectrum, as a reference in order to check for weak, isolated lines. 

Results. We test the quality of the new line list by re-deriving stellar parameters for 45 1 stars with high resolution and signal-to- 
noise HARPS spectra, that were analyzed in a previous work with a larger line list. The comparison in temperatures between this 
work and the latest IRFM for the stars in common shows that the differences for the cooler stars are significantly smaller and more 
homogeneously distributed than in previous studies for stars with temperatures below 5000 K. Moreover, a comparison is presented 
between interferometric temperatures with our results that shows good agreement, even though the sample is small and the errors of 
the mean differences are large. We use the new line list to re-derive parameters for some of the cooler stars that host planets. Finally, 
we present the impact of the new temperatures on the [Cr i/Cr n] and [Ti i/Ti n] abundance ratios that previously showed systematic 
trends with temperature. We show that the slopes of these trends for the cooler stars become drastically smaller. 

Key words, techniques: spectroscopic - stars: fundamental parameters 
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1. Introduction 

Temperature (T e g), surface gravity (logg), and metallicitity 
([Fe/H], where iron is used as a proxy) are basic atmospheric 
parameters necessary to characterize a star, as well as to deter- 
mine other indirect and fundamental parameters, such as mass, 
radius and age that are acquired in combination with stellar evo- 
lutionary models ( |Girardi et al.|2000) . Precise and accurate stel- 
lar parameters are also essential in exoplanet searches. The light 
curve of a transiting planet orbiting a star gives information of 
the planetary radius always in dependence of the stellar radius 
(R^ oc R + ). Moreover, the mass of the planet, or the minimum 
mass in case the inclination of the orbit is not known, is de- 
termined from the radial velocity technique only if the mass of 
the star is known (M p oc MJ ). Therefore, the determination of 
the planetary radius and mass requires a combination of tran- 



siting and radial velocity data (Ammler-von EifF et al. 2009 
Torres et al. 2008, 2012). Apart from the derivation of planetary 
properties (mass, radius, density), stellar parameters can be used 
to reveal correlations between planets and their host stars that 



* Based on observations collected at the La Silla Paranal Observatory, 
ESO (Chile) with the HARPS spectrograph at the 3.6-m telescope (ESO 
runs ID 072.C-0488. 



will give insights on their formation and evolution mechanisms. 
Several correlations have been studied so far with very inter- 
esting results, such as stellar metallicitity and planet frequency 
(e.g. |Santos et al.|2004|[Udry et al.|2007t|Sousa et al.|201 1) , stel- 
lar metallicitity and planetary mass (Guillo fet al.|200 6), stellar 
metallicitity and planetary orbital periods (|S ozzetti 2004), stel- 
lar temperature and obliquities ( Albrecht et al. 2012; Winn et al. 
2010), metallicitity and planet radius (Buchha ve et al.||2012) . 
It is obvious that in order to correctly characterize the planets 
and furthermore, to statistically verify such correlations, well- 
determined stellar parameters are required. 



Hundreds of stellar spectra are available from radial veloc- 
ity planet search programs (Udry et al. 20 00} Wgt et al.||2 000 



|Mayor et al.pOOlj [Bouchy et al.||2009[ |Lo Curto et al.||2010) . 
These spectra are obtained with high resolution spectrographs 
and a combination of them gives high signal-to-noise (S/N) that 
makes the spectroscopic analysis more powerful. 



The study of the Ha wings (Fuhrmann 200 4JI, the e xcitation 
and ionization balanc e of iron lines (jSantose t al. 2004), spectral 
synthesis techniques ( Valenti & Fischer 2005), line ratios (or line 
depth ratios) (Gray 1994) are the basic spectroscopic techniques 
that along with fundamental techniques such as photometry and 
interferometry can be used for the determination of the effective 
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temperature. A comparison between these different methods can 
show considerable discrepancies in their results. Even in the re- 
stricted group of solar-type stars, the effective temperatures ob- 
tained with these methods can differ significantly (e.g.|Kovtyukh| 



Table 1. Characteristics of the sample and the reference star, 
described in SO08. 



et al.|2003[ [Ramirez & Melendez|2004"t|Casagrande et al.|2006 
Sousa et al.|2008| >. 

The temperature determination becomes more difficult when 
we focus on K-type stars. The difficulties in these stars with 
T e ff < 5000 K emerge from their line crowed spectra that cause 
strong blending. Blending can be a considerable problem if one 
uses the standard technique based on the iron equivalent widths 
(EWs). Lines cannot be easily resolved and the continuum place- 
ment becomes more difficult, causing bad measurement of the 
EWs and hence, makes the calculation of stellar parameters am- 
biguous. Therefore, it is important to select carefully the iron 
lines in such manner that will eliminate the blending effects, es- 
pecially for cool stars. In addition, the choice of the atomic pa- 
rameters influences the abundance determination. Some authors 
calculate the atomic parameters using the Sun as a reference to 
avoid the errors that emerge from the theoretical or laboratory 
values. In that way, the atomic parameters for stars that are dif- 
ferent from the Sun, i.e. too hot or too cool can be no longer 
accurate enough. 

The accuracy in the fundamental parameters of planet host 
stars is of great importance to planetary studies since these er- 



rors can propagate to the planetary properties (e.g. Pepe et al. 



201 1 1. In addition, stellar parameters affect the determination of 
other element abundances ([X/H]) and can introduce potential 
biases in their abundance calculations (e.g. Neves et al.| J2009 
Adibekyan et al.|2012b| > 



In this paper, we use the data from previous work that per- 
formed the standard spectroscopic analysis for a sample of solar- 



type stars presented in the paper of Sousa et al. ( 2008 ) (hereafter 
SO08). This sample is part of the High Accuracy Radial ve- 
locity Planet Searcher (HARPS) Guaranteed Time Observations 
(GTO) survey that is composed of slow rotators and low activity 
FGK stars in order to detect low mass planets. 

A comparison of these spectroscopic results with the Infrared 
Flux Method (IRFM) indicates a disagreement between the ef- 
fective temperatures only for the cooler stars of the sample with 
temperatures below ~ 5000 K. Motivated by that, we compile an 
optimized line list to improve the accuracy of the stellar param- 
eters for the cooler stars and compare our results with other in- 
dependent methods (IRFM, interferometry). This paper is orga- 
nized as follows: in Sect. [2] we describe the stellar sample, the 
spectroscopic method and the criteria to select an optimized line 
list. We compare our results with the previous work of SO08 in 
Sect. [3] In Sect. HI we compare the spectroscopic surface gravi- 
ties with the trigonometric ones and address the ionization prob- 
lem for the cooler stars of our sample. In Sect.B] we compare the 
new derived parameters with the IRFM and interferometry. In 
Sect. [6] we calculate new parameters for some cool planet hosts. 
Finally, in Sect. [7] we use the new parameters for the cooler stars 
to reduce the trends of other elements ([Cri/Cru] and [Tii/Tin] 
abundance ratios) with temperature, as reported in Adib ekyan] 
|etaLl ( 12012bD . 



2. Stellar sample and previous spectroscopic 
analysis 

The stellar sample, presented in SO08, is composed of 451 stars 
as part of the HARPS high-precision GTO program at the ESO 
La Silla 3.6m telescope with the objective to detect low-mass 



Lowest 
Highest 
Reference star 



r eff logg [Fe/H] Mass 

(AT) (dex) (dex) M 

4556 3.68 =084 037" 

6403 4.62 0.39 1.42 

4723 4.40 -0.02 0.76 



extra-solar planets with high radial velocity accuracy (Mayor 
et al. |2003] l . It is mainly comprised of dwarf FGK stars selected 
from a volume-limited sample of the CORALIE survey ( Udry 
et al. 2000). Planet hosts from the southern hemisphere were 
also added to this sample, forming in total a sample of 451 stars. 
These stars are slowly-rotating, non-evolved, and low-activity 
stars, with apparent magnitudes that range from 3.5 to 10.2 and 
have distances of less than 56 parsec. The spectra have a resolu- 
tion of R ~ 1 10,000 and 90% of the combined spectra have S/N 
higher than 200. We point to SO08 for more details. 

For this sample, SO08 derived stellar parameters by impos- 
ing excitation and ionization equilibrium, based on the measure- 
ments of weak iron lines. This method is very effective for FGK 
stars due to the numerous iron lines in their spectra. Iron abun- 
dance is used as a proxy for the overall stellar metallicitity. 

The line list for their spectroscopic analysis, was composed 
of 263 Fei and 36 Fen lines. The EWs of the lines were 
measured automatically for all stars using the ARESr] code 
(Automatic Routine for line Equivalent widths in stellar Spectra; 
Sousa et al.|[20 07). The atomic parameters of the iron lines, 
namely the oscillator strength values (log gf), were computed 
by an inverted solar analysis, using a solar model with T e g - 
5777 K, logg = 4.44 dex, &= 1.0 m s~\ \og e (Fe) = 7.47 dex. 

The spectroscopic analysis was completed assuming Local 
Thermodynamic Equilibrium (LTE), and using the 2002 version 



(Sneden 1973) 



of the abundance determination code MOOGr 
and a grid of Kurucz Atlas 9 plane-parallel, ID static model at- 
mospheres (Kuruc zjl993| ). 

The best parameters were obtained when the Fe i abundance 
shows no dependence on the excitation potential, xu and on 
the reduced equivalent width, log W^/A. Additionally, the mean 
abundances given by Fei and the Fen must be the same (ion- 
ization balance) and consistent with those of the input model at- 
mosphere. The input parameters converge to the true ones with 
an iterative minimization code based on the Downhill Simplex 
Method (Press et al.||1992|), making the total procedure auto- 



matic ( |Santos et al. 2004 1. Some characteristics of the sample 
are depicted in Table 111 as described in SO08. 

There are different sources of uncertainties that occur in the 
stellar parameter determination using this method. These errors 
can be attributed to the uncertainties of the measurements of the 
EWs, the uncertainties in the atomic parameters and the uncer- 
tainties that are intrinsic to the method of ionization and excita- 
tion equilibrium. In addition, systematic errors can arise due to 
the assumptions of the method, such as ID static atmospheres, 
NLTE effects ( |Mashonkina et al.|201 l] |Bergemann et al.|2012) . 
However, departures from LTE for Fe lines do not need to be 



1 The ARES code is an open source code and can be found at 
http://www.astro.up.pt/~sousasag/ares 

2 The MOOG code can be downloaded free at 
http ://verdi .as.utexas. edu/moog .html 
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considered for near solar metallicity dwarfs but should be taken 
into consideration for more evolved or very metal-poor stars 
(Lind et al. 2012 Ruchti et al. 2013 1 that are not part of this 
sample. 

Errors in the measurements of the EWs can be minimized 
by using high quality spectra. In low S/N spectra, weak lines 
cannot be distinguished from noise and strong lines can be un- 
derestimated due to the miscalculation of their wings. The high 
resolution and high S/N spectra used for this sample, are the best 
solution to deal with such errors. Since in our spectroscopic anal- 
ysis the atomic data (log gf) are derived with respect to the Sun, 
we expect small errors for solar analogs and more significant for 
cooler and hotter stars. 

The errors in the stellar parameters for this analysis are cal- 
culated by varying each parameter (temperature, surface gravity 
and microturbulence) by a standard value. The uncertainty in £, 
is determined from the standard deviation in the slope of the 
least-squares fit of log e (Fei) versus log Wa/ A. The uncertainty 
in T e ff is determined from the uncertainty in the slope of the 
least-squares fit of log e (Fe i) versus xi m addition to the uncer- 
tainty in the slope due to the uncertainty in £,. The uncertainty in 
log g is derived from the contribution of the uncertainty in T e g 
and microturbulence in addition to the scatter error of the Fe n 
abundance (measured as <x/ Vw, cr is the standard deviation and 
N the number of lines). The uncertainty in the log e (Fe i) is the 
sum of the squared uncertainties due to the error in r e ff , f f , and 
the scatter error of Fe i abundance. The use of many iron lines 
can reduce this type of uncertainty, assuming that the majority 
of the lines are independent and of good quality. 



2. 1 . Building a stable line list for the cooler stars 

A reliable line list is comprised of lines that can be accurately 
measured, which usually means unblended lines. In addition, 
lines must be unsaturated, cover a wide range in excitation po- 
tential and have accurate atomic data. Temperature, as well as the 
other stellar parameters, is strongly correlated with the equiv- 
alent width. This sensitivity emerges from the excitation and 
ionization processes that follow the exponential and power de- 
pendencies with temperature that are defined by the well-known 
Boltzmann and Saha equations. 

For the cooler stars line blending is severe, which makes the 
measurements of the EWs problematic. In particular, blending 
effects cause an overestimation of the EWs as two blended lines 
cannot be resolved. 

Another bias in the EW measurements may come from the 
fitting of strong lines. Gaussian fitting is a good approximation 
for weak lines and it can be reliable up to 200mA based on our 
experience, whereas a Voigt profile should be used for stronger 
lines. Saturated lines that deviate significantly from the linear 
part of the curve of growth should also be avoided in the abun- 
dance analysis. The EW predicted by the models of strong lines 
that are highly saturated, is quite dependent on microturbulence. 
A wrong estimation of microturbulence, will then produce errors 
in the abundance of any highly saturated line. 

On the other hand, weak lines that are strongly blended could 
lead to a sub-estimation of the continuum and consequently of 
the EW. This effect, however, is less significant. An overestima- 
tion in the EWs due to blending, as well as the underestimation 
of very strong lines could be the reason for the systematic raise 
in temperature that is observed for the cooler stars of SO08. In 
addition, the reduced equivalent width could also be affected by 
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Fig. 1. Upper panel: Curve of growth for both line lists for the ref- 
erence star computed for the temperature (0=5O4O/r eff ) of this work. 
Circles represent the reduced EW of the line list of this work and crosses 
the line list of SO08. Lower panel: The Fe i abundances of the reference 
star versus the excitation potential. The dashed line shows the positive 
slope that corresponds to the line list of SO08. The solid line corre- 
sponds to the line list of this work and the slope is obviously zero. 



Table 2. Sample of the line list used for the spectroscopic anal- 
ysis with the atomic parameters of Fe i and Fe n as well as the 
corresponding EWs of the reference star HD 21749. 



A (A) 


Xi 


logg/ 


Element 


EW(mA) 


4508.28 


2.86 


-2.403 


Fen 


53.0 


4520.22 


2.81 


-2.563 


Fen 


72.4 


4523.40 


3.65 


-1.871 


Fel 


101.7 


4537.67 


3.27 


-2.870 


Fel 


43.2 


4551.65 


3.94 


-1.928 


Fel 


41.9 


4556.93 


3.25 


-2.644 


Fel 


57.9 


4566.52 


3.30 


-2.156 


Fel 


68.6 


4574.22 


3.21 


-2.353 


Fel 


55.1 


4576.34 


2.84 


-2.947 


Fen 


29.6 



such biases, leading to correlations with the excitation potential 
(see Appendix \A\. 

Therefore, our aim is to optimize the iron line list of SO08. 
With this goal, we use the K-type dwarf HD 21749 with T e s 
= 4723 K (see Table [TJ, as reference in order to check for un- 
blended lines in its high S/N spectrum. After visual inspection, 
we only consider weak, isolated lines that can give good esti- 
mation for the local continuum. We avoid strong lines (>200 
mA) in order to apply Gaussian profiles. For the reference star 
we show the curve of growth (see |Gray|[T994 1 using both line 
lists (Fig.[T]upper panel). Limiting the EW cut off, we mitigate 
in large amount the problem of saturated lines and microturbu- 
lence. The proof of that mitigation is the fact that the derived 
temperatures with the new line list agree with other less model 
dependent methods (see Sect.Bl). 

Very weak lines (<10 mA) were also excluded so that noise 
is not superposed to these lines. The region of the spectrum be- 
low 4500 A is neglected due to the higher blending. The final 



M. Tsantaki et al.: Deriving precise parameters for cool solar-type stars 



Table 4. Results of the internal comparison for the whole sample 
and for stars with r e ff<5000K. 



Ar cff 

(K) 



Alogg 
(dex) 



A [Fe/H] 
(dex) 



This Work - SO08 
whole sample 

This Work - SO08 
T eff <5000K 



-31+3 -0.023±0.003 0.0002+0.0011 
-106±6 -0.07+0.01 0.0127+0.0001 



line lisfjis compiled with 120 Fe i and 17 Fe n lines, as shown in 
Table ET 

As mentioned before, the effective temperature is derived 
when the correlation coefficient between log e (Fe i) and xi is 
zero. In the lower panel of Fig. [T] we demonstrate this corre- 
lation for the reference star using the line list of this work and of 
the work of SO08 using the parameters derived with the line list 
of this work. The positive slope for the line list of SO08 is trans- 
lated in an overestimation in temperature of ~180 K for this star. 
In addition, the abundances with the new line list show a smaller 
scatter which corresponds to smaller errors in the final tempera- 
ture value. 



3. New stellar parameters for 451 FGK stars in the 
HARPS GTO sample 

In order to check the effectiveness of the new line list, we re- 
derive stellar parameters for the 451 stars of the sample. For 
consistency, we use the same EWs as in SO08 that were mea- 
sured automatically for all stars with the ARES code. In addi- 
tion, we use the same damping parameters and atomic data. After 
a preliminary determination of the fundamental parameters, we 
perform a '3cr clipping' procedure for lines that contribute with 
abundances higher than 3cr from the average abundance. This 
procedure was also applied in SO08. 

TablePJlshows the final parameters for a fraction of stars from 
the sample. Microturbulence (£,) is used as a free parameter and 
is also derived from this spectroscopic analysis. The correlation 
of microturbulence with temperature and surface gravity is pre- 
sented in Appendix [B] This calibration can be useful in cases 
where the value of £, is set fixed. The stellar masses are cal- 
culated using the stellar evolutionary models from the Padova 
grourjj The errors of the fundamental parameters are internal, 
attributed to the method. They, thus, represent relative errors and 
not the absolute accuracy. 

3. 1 . Internal comparison 

We present the comparison results with the work of SO08 for 
temperature, surface gravity and metallicitity, respectively in 
Fig. [2] Temperatures show very good agreement in the high 
and intermediate ranges. Table [4] shows the mean difference in 
Jeff for the whole sample and for temperatures below 5000 K 
with their standard errorqj The mean difference in r er y is -31+3 
(cr=53) K for the whole temperature range. The significant dif- 
ferences appear, as expected, for stars with temperatures below 
5000 K with mean difference Ar eff = -106±6 (<r=54) K. 



3 The full line list is available in electronic form. 

4 Web interface for stellar mass estimation: 
http://stev.oapd.inaf.it/cgi-bin/param 

5 The standard errors of the mean (<x M ) are calculated with the fol- 
lowing formula: cr M =-j=, <x being the standard deviation. 
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Fig. 2. Comparison between the parameters derived with the cool line 
list of this work and the results of SO08 for: temperature (top panel), 
metallicitity (middle panel) and surface gravity (bottom panel). Ar cff 
corresponds to this work minus SO08. Black squares represent stars 
with r cff <5000 K. Triangles represent stars with planets taken from 
Table|7](see Sect.|6j. 



One interesting result is that even though we have 
considerably large differences in the low temperature 
regime, the values of metallicitity remain unaffected with 
< A[Fe/H] >=0.00+0.00 (cr=0.02) dex for the overall sample. 
The same effect appears for surface gravity with < Alogg >=- 
0.02+0.00 (cr=0.06) dex, even though there is bigger scatter. 
This result suggests that surface gravity and metallicity are not 
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5^ 









Star 


T e g 


logs ?< 


[Fe/H] 


Mass 


Age 








(K) 


(cm s~ 2 ) (kms -1 ) 


(dex) 


(Mo) 


(Gyr) 




HD967 


5595+17 


4.59+0.02 0.90+0.05 


-0.66+0.01 


3.76+0.02 


7.8+3.9 




HD 1237 


5489+39 


4.46+0.11 1.04+0.06 


0.06+0.03 


3.92+0.02 


1.8+1.7 




HD 1320 


5699±12 


4.55+0.05 0.89+0.02 


-0.26+0.01 


3.91+0.02 


2.4+2.0 




HD 1388 


5970±15 


4.42+0.05 1.13+0.02 


0.00+0.01 


1.06+0.01 


3.0+0.8 




HD 1461 


5740±16 


4.36+0.03 0.94+0.02 


0.18+0.01 


1.05+0.02 


2.0+1.1 




HD 1581 


5990+15 


4.49+0.07 1.24+0.03 


-0.18+0.01 


1.03+0.02 


2.2+1.0 




HD 2025 


4851+49 


4.49+0.13 0.51+0.18 


-0.37+0.02 


3.71+0.01 


4.6+4.0 




HD2071 


5729+12 


4.49+0.02 0.93+0.02 


-0.08+0.01 


3.99+0.01 


0.9+0.6 




HD 2638 


5169+77 


4.41+0.15 0.66+0.13 


0.12+0.05 


3.85+0.03 


4.0+3.8 
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Fig. 3. The effect of temperature on the other parameters: metallicitity (left panel) and surface gravity (right panel). The differences correspond to 
the values of this work minus the results of SO08. 



as sensitive to the selection of the line list as temperature, for 
this temperature regime and for this method. The same effect 
appears for the cool stars as seen in Table HI 



The impact of the updated effective temperatures on [Fe/H] 
and on \ogg is depicted on Fig. [3] The changes in metallicitity 
show almost no correlation, within the errors, with the effective 
temperature. Only slightly higher metallicities appear for low 
temperatures, when comparing with SO08, yet within the errors. 
The surface gravities are also not correlated with temperature, 
even though there is high dispersion in the low temperature re- 
gion. This result suggests that using this technique, a potential 
error in one of the parameters will not propagate to the others, 
avoiding systematic errors (see also [Torres et al.|2012| ). 



Metallicity has a key role in planet formation theories and is 
correlated with the planet frequency. The stellar sample, as men- 
tioned before, is part of the HARPS GTO planet search program 
and it contains 102 up-to-date planet hosts. The metallicitity dis- 
tribution of the sample, presented in SO08, shows that the Jovian 
planets are preferentially found in metal-rich stars, in contrast to 
Neptune-like planets that do not seem to follow this trend, even 
though the number of these planets is small. The new metal- 
licities derived with the new line list do not change this trend, 
making this correlation between stars and planets reliable even 
before adapting the new temperatures. 



4. Ionization balance problem in cool stars. 

In an LTE abundance analysis for solar type stars the ioniza- 
tion equilibrium should be satisfied. Using the standard spectro- 
scopic method we force the abundance of Fe i and Fe n to agree. 
Surface gravity is determined from this tuning. However, there 
are many studies for cluster stars ( Yong et al.||2004 Morel & 
Micela 2004; Schu ler et al |2006]l and fi eld stars ( |Allende Prieto 
et al. 2004; Ramirez et al. 2007, 2013) showing that the cooler 
dwarf stars deviate from the ionization balance, with system- 
atic higher Fe n abundances over Fe i. The authors explain these 
discrepancies due to possible different scales in the stellar pa- 
rameters (namely T e s and log g) or due to NLTE effects caused 
by the simplifications of the model atmospheres. The ionization 
balance of Fe i and Fe n can be investigated in this work by the 
behaviour of surface gravity. In fact, surface gravity mostly de- 
pends on Fe ii lines, once the temperature scale is correct. 

An essential test for the accuracy of surface gravity derived 
from spectroscopy is the comparison with surface gravity de- 
rived from parallaxes (trigonometric log g), based on the funda- 
mental relation: 

g M T e ff L 

log ±- = log — + 41og —^- - log — 

80 % T eff,& L 

We calculated the trigonometric \ogg using the new 
Hipparcos parallaxes ( |van Lee uwen 2007 ), V magnitudes, bolo- 
metric correction based on Flower (1996) and |Torres| ( |2010| l, 



solar magnitudes from (Bess ell et al.||1998) , the spectroscopic 
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Fig. 4. Comparison of the surface gravities derived from spectroscopy, 
log g spec of this work and from Hipparcos parallaxes, logg BIP . A mean 
error for both axes is given in the upper left part. 
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Fig. 5. Trigonometric log g HIP minus logg spec as a function of temper- 
ature. The dashed line represents a linear fit (-2.578T0 _4 r cf f+1.477). 
The mean error is shown at the bottom left. 



masses and r e ff. No correction for interstellar reddening is 
needed since all stars are less than 56 pc in distance. 

In Fig. HI we compare the spectroscopic log g spec with the 
trigonometric log gnip- At first glance, the spectroscopic logg 
agrees with the trigonometric (< \oggnip - \ogg spe c >=0.07 
dex). However, a more careful look shows that there is a dis- 
agreement especially for the high values of logg (>4.5 dex) 
where the spectroscopic gravities are underestimated, which is 
also observed in the work of SO08. 

An interesting fact is that the differences between the 
trigonometric and spectroscopic log g are greater for the cooler 
stars (Fig. Bl). There is a clear trend between the differences in 
log g and r e ff , where the underestimation of surface gravity in 
low temperatures becomes higher (< log gHip - log g spec >=0.22 
dex for stars with r er r < 5000 K). This is translated into system- 
atically higher Fe n abundances over Fe i for low T^, 

Such differences between Fe i and Fe n abundances are dif- 
ficult to explain with model uncertainties and departures from 
LTE in the spectral line formation calculations as they are ex- 
pected for the warmer stars of our sample where most iron is 
ionized ( |Lind et al.|2012[ ). In addition, other model uncertainties 
related to granulation and activity of K-stars have been proposed 



Fig. 6. Comparison between the temperatures derived from this work 
and the IRFM for stars in common. 



to explain these differences even though these effects should be 
evident for young stars (Morel & Micela 2004 Schuler et al. 
[20T0| . 

Other possible explanations for these differences have to do 
with the iron ionization method itself. We use Fen lines that 
strongly depend on surface gravity. For solar-type stars, how- 
ever, these lines are not sufficiently present leading to poorly 
constrains of log g. On the other hand, the numerous Fe i lines 
are insensitive to logg changes. 

We have to note though, that temperatures and metallicities 
derived using the ionization and excitation equilibrium of iron 
lines are shown to be mostly independent of the adopted surface 
gravity ( [Torres et al.|2012 i. Hence, the temperatures and metal- 
licities derived with our spectroscopic method can be used as 
reference even if the derived spectroscopic surface gravities dif- 
fer from the trigonometric values (see also Santos et al.| 20137 



subm itted) . 



5. Comparison with other methods 

In order to evaluate the consistency of our results, namely for 
r e ff, we compare them with other techniques. Here, we present 
a comparison with two different methods that are considered to 
be less model dependent, the Infrared Flux Method (IRFM) and 
interferometry, respectively. 



5. 1 . The Infrared Flux Method - IRFM 



The IRFM method ( |Blackwell & Shallis|197"7) is a semi-direct 
method for determining stellar parameters. The principle of this 
method relies on the fact that the bolometric flux depends on the 
angular diameter and the effective temperature, as described by 
the Stefan-Boltzmann law, whereas the monochromatic flux in 
the infrared (IR) depends on the angular diameter but weakly on 
the effective temperature, this way the dependence on the angu- 
lar diameter disappears: 



fl 



hoi 



o-T 4 



I air f,\iR{model) 



(1) 



where //,„/ is the measured bolometric flux, /xir is the mea- 
sured monochromatic IR flux and f /UR (model) is the monochro- 
matic flux in the IR derived by the assuming model. The IRFM 
has the advantage that the dependence on the models is limited 
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Table 5. Comparison between the effective temperatures derived 
with different methods, <x represents the standard deviation and 
N the number of stars for the comparison. 



Method 


AT ef[ (K) 


<r(K) 


N 


IRFM - this work 
IRFM - SO08 


+33 ±3 
+ 14±3 


54 
61 


347 
347 


IRFM - this work (T cff <5000K) 
IRFM - SO08 (r eff <5000 K) 


+42 + 12 
-86+16 


65 

86 


29 
29 


Interferometry - this work 
Interferometry - SO08 
Interferometry - IRFM 


-4 ±33 
-56 ± 35 
-98 ± 83 


98 

105 
204 


9 
9 
6 



while the spectroscopic effective temperatures have considerable 
model dependence. 



We compare our results with the work of Casagrande et al 



(2010| |2011) that implement the IRFM method for a large sam- 
ple of stars. The authors estimate the bolometric flux from multi- 
band photometric measurements in the optical BV(RC)c band 
and in the near-IR 2MASS JHK$ band. For the missing spec- 
tral regions, the flux is calculated by synthetic spectra computed 
from model atmospheres. The absolute calibration of Vega is 
based on its synthetic spectrum with an uncertainty of the zero 
point of -15 K ( |Casagrande et al.|2010) . 

Figure [6] depicts the comparison between the spectro- 
scopic temperatures and the IRFM for the stars in common. 



Temperatures of 341 stars are taken from Casagrande et al 
(201 1 ) using stars with direct application of the IRFM (irfm sam- 
ple) and stars with T e ff derived from colour calibrations (clbr 
sample). Moreover, 6 stars were taken from Casagrande et al. 
(2010). The comparison between the results of this work and the 
IRFM shows good agreement for all temperature ranges. In par- 
ticular, for the cooler temperature region, the differences in T e g 
between this work and the IRFM are much smaller and more 
homogeneously distributed than between SO08 and the IRFM. 

The mean differences in temperature for the comparison 
samples are shown in Table [5] It is clear that the differences in 
temperature for this work with the IRFM are constant through- 
out the temperature range, with a small offset of 33 K for the 
whole sample. For the cooler stars these differences are AT e s = 
42 ± 12 K that are much smaller than SO08 with AT eS = -86 ± 
16 K. Figure 17] shows the comparison of stars with r e ff<5000K. 
It is evident that the trend in Ar e ff of SO08 mostly disappears 
with the new temperatures. 

5.2. Interferometry 

Precise measurements of stellar angular diameters are acquired 
through long baseline interferometry. The standard practice to 
determine the angular diameter is to fit the observed visibilities 
as a function of baseline to a uniform disk model. The angular 
size is connected to the more realistic limb darkened angular size 
(Old) using correction factors from model atmospheres (Claret 
2000). Temperature is then derived with the standard relation: 
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Fig. 7. Comparison between the difference in temperatures derived 
from IRFM - This Work (circles) and IRFM - SO08 (triangles). The 
dashed and the solid line depict the linear fits of the data with slopes: 
+0.20+0.04 and +0.43+0.05, respectively. 
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where fs i is usually calculated from the Spectral Energy 
Distribution. 



Fig. 8. Comparison between the spectroscopic, IRFM and the direct 
temperature measurements for stars in common with the literature. The 
x-axis corresponds to temperatures derived from this work. Square sym- 
bols are for the different temperatures for HD 26965, stars for HD 10700 
and triangles for HD 146233. Stars with grey colour represent T e g with 
high angular diameter uncertainty. 



We compare our results with the temperatures derived from 
interferometry. Unfortunately, the number of stars with avail- 
able angular diameters for this sample is only down to a few 
since these measurements are challenging for dwarfs due to their 
small photospheric disks that are difficult to resolve. We have 9 
stars in common for the comparison with spectroscopy and 6 
with the IRFM. We use only direct angular diameters and bolo- 
metric fluxes available in the literature from Table [6] Our re- 
sults show better agreement for these stars to interferometry than 
when comparing with the values of SO08 and the IRFM (see 
Table [5]). We have to note though, that the comparison sample 
is very small and the values of r e ff were derived from the clbr 
sample that is not the best representative for the IRFM precision. 

The differences in temperatures for the different methods are 
plotted in Fig. [8] For 3 stars (HD 10700, HD 26965, HD 146233) 
we include angular diameter measurements from different au- 
thors that give different r e ff and are represented with different 
symbols. In the same figure, we see that from the stars with 
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Table 6. Interferometric data and derived temperatures for stars in common with our sample. 



Star 
HD 



10700 

11964 
19994 
22049 
23249 
26965 

128621 
146233 

209100 



8ld 
(mas) 



A6/6 
(%) 



1 elT 

(K) 



(K) 



r rtlus\ 
1 eff 



(K) 



' eff 

(K) 



References 



2.022 : 

1.971 : 
0.611: 

0.788 : 

2.148: 

2.394 : 

1.504: 
1.650: 
6.001 : 

0.676 : 

0.780 : 
1.890: 



0.011 
0.050 
0.081 
0.026 
0.029 
0.029 
0.006 
0.060 
0.021 
0.006 
0.017 
0.020 



0.54 
2.54 
13.25 
3.30 
1.35 
1.21 
0.40 
3.63 
0.35 
0.89 
2.18 
1.06 



5383 ± 47 
5449 ± 83 
5413 ± 359 
6109 ± 111 
5107 ±21 
4986 ± 57 
5143 ± 14 
4910 ± 90 
5182 ±24 
5836 ± 46 
5433 ± 69 
4527 ± 29 



5310 ±17 5322 ±17 5459 ± 80 



5332 ± 22 
6289 ± 46 
5153 ±42 
5150 ± 51 
5153 ±38 



5285 ± 21 
6315 ±44 
5049 ± 48 
5027 ± 48 
5098 ± 32 



6159: 

5207: 



80 
80 



5311 ±80 



5234 ±63 5168 ±75 

5818 ±13 5810 ±12 5826 ± 80 

4754 ±89 4649 ±73 4731 ± 80 



l,a 

2, a 
3,b 
3,b 

4, c 

5, a 
6,d 
4,d 

7. e 

8. f 

9. f 
4,g 



No tes. Referenc es for 9 LD : (1) Teixeira et al. (2009); (2) Pijpers et al. (20 03}; (3 ) [van Belle & vonBr aun (2009); (4) Kervella & Fouque (2008) 



(5) Thevenin et al. (2005 1; (6) Boyajian et al. (2012b); (7) Kervella et al. (2003); (8) 
for /fa, ,: (a)|Bruntt et al. (2010); (b) van Bell e & von Bra un ( 2009 l;~(c) lCayrel et al.| ( [20TT 
pOOl) ; (f) |Boyajian et al.| ( |2012a) ; (g) |Ramirez & Melendez|j2005[ ). 



Bazot et al.H2011}; (9 ) Boyajian et al. (2012a). Refer ences 
(d) |Boyajian et al.| < |2012b) ; (e) |Ramirez & Melendez] 



multiple measurements, the ones with the smallest uncertainty 
(^%) agree better with the temperatures derived with the spec- 
troscopic and photometric methods. A precision better than 2% 
in the angular diameters corresponds to an accuracy of 1 % in the 
effective temperatures, which is roughly 60 K at solar tempera- 
ture, assuming no error in the bolometric flux. It is useful thus, 
to take the uncertainty of the angular diameter into consideration 
for a reliable determination of temperature. 



6. New atmospheric parameters for cool planet 
hosts 

In general, for planet host stars with effective temperature below 
5200 K, the new set of parameters, derived using the line list pre- 
sented in this paper, imply that the planet host stars have lower 
temperatures than previously published. This has implications 
for both their mass and radius determination. The lower temper- 
atures imply lower stellar masses as well as lower stellar radii. 
This means that the derived planetary masses and radii (for tran- 
sit planet cases) are also lower. As a consequence of the mass 
reduction, the semi-major axis of the orbits will also be smaller. 
The expected effects are however small, and no major revisions 
are expected to occur. 

With this new line list, we re-derive the stellar parameters for 
10 "cool" planet hosts already published in the literature and are 
not included in the 45 1 stellar sample of this work. We only con- 
sider GK dwarfs with an effective temperature lower than 5200 K 
whose planets were detected with the radial velocity technique 
from the CORALIE and HARPS GTO planet search samples. 
These planet hosts have been previously analyzed with high S/N 
spectra following the same procedure as this work but with dif- 
ferent line lists. 

In TablelT] the fundamental parameters based on the new line 
list are presented. The sixth column gives the reference of the 
previously published parameters. To explain the effect of these 
new parameters to mass more quantitative, we calculate the stel- 
lar masses from the Padova interface for both with the origi- 
nal and new parameters. We avoid to use the published stellar 
masses in order to compare uniformly. We find the maximum 
difference in mass to be 1.5% in absolute units which is negligi- 
ble compared to the standard mass error. 



7. [Cr i/Cr n] and [Ti i/Ti n] vs. r e ff with the new 
atmospheric parameters 

The precise and accurate stellar parameters are, as well, very im- 
portant for further analyzing stellar chemical abundances. The 
traditional spectroscopic abundance analysis methods require 
these parameters as input to compute the atmosphere models, 
hence the accuracy of the final elemental abundances depends 
on the accuracy of these input parameters. Different atoms and 
ions are not equally sensitive to all the stellar parameters. For 
example, ionized species are more sensitive to gravity variations 



2006 Neves et al. 2009 



than neutral species (e.g. Gilli et al 
Adibe kyanet al |2012b). 

Recently, |Neves et al.| ([2009 ) and Adibeky aiTet al.| ( |2012b] > 
analyzing chemical abundances of the refractory elements of the 
HARPS sample stars, observed some unexpected trends with ef- 
fective temperature. Particularly, they detected systematic trends 
of [X/H] or [X/Fe] with r e ff for some elements at low temper- 
atures and found that [Cr i/Cr n] and [Ti i/Ti n] abundance ratios 
gradually increase with decreasing effective temperature when 
r e ff < 5000 K. Similar trends for different elements with T e g 



have been already noted in the literature (see e.g. Valenti & 
Fischer 2005; Pres ton et al.||2006[ IGilli et~dl|2006| |Lai et al. 
2008 ; Suda et al. 2011). Different explanations of the mentioned 
trends are discussed in the literature. The unexpected trends in 
the low temperature regime may be due to the stronger line 
blending and may also be connected to either deviations from ex- 
citation or ionization equilibrium, or to problems associated with 
the differential analysis ( [Neves et al.|20 09). A possible explana- 
tion for the observed trends with T e g could also be an incorrect 



T-t relationship in the adopted model atmospheres (Lai et al. 
2008) or NLTE effects ( [Bodaghee et al.|2003) >. Summarizing, it 
can be assumed that the observed trends are probably not an ef- 
fect of stellar evolution, and uncertainties in atmospheric models 
are the dominant effect in measurements (see also the discussion 
in |Adibekyan et al.|20i2"b"| ). 

In Fig. [9] we plot the [Cri/Cru] and [Tii/Tin] abundance 
ratios derived using the stellar parameters of SO08 and this 
work as a function of the r e ff for stars cooler than 5000 K. 
This plot is useful to ensure that the ionization equilibrium en- 
forced on the Fe n lines is acceptable to other elements. As can 
be seen the slopes of the abundance ratios with new parame- 
ters are very gentle. The new slope of [Cri/Cru] per 1000 K is 
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Table 7. Updated stellar parameters for previously analyzed planet hosts. 



Name 



T e ff 

(K) 



lOj 

(cm 



',g 

s- 2 ) 



[Fe/H] 
(dex) 



(km s" 1 ) 



BD-082823 
HD3651 
HD 13445 
HD 20868 
HD 99492 
HD 125595 
HD 128311 
HD 192263 
HD 215497 
HIP 5158 



4648 ± 135 
5182 ±79 
5114 ± 61 

4720 ± 91 
4815 ±184 
4596 ± 235 
4778 ± 75 
4906 ± 57 
5003 ± 103 
4673 ± 175 



4.33: 
4.30: 
4.55: 
4.24: 
4.28: 
4.25: 
4.35: 
4.36: 
4.26: 
4.24: 



0.32 
0.16 
0.13 
0.22 
0.46 
0.63 
0.17 
0.17 
0.26 
0.47 



0.00 ± 0.08 
0.12 ±0.05 
-0.29 ± 0.04 
0.08 ± 0.06 
0.24 ±0.12 
0.10 ±0.14 
-0.03 ± 0.02 
-0.07 ± 0.02 
0.25 ± 0.05 
0.22 ±0.12 



0.27: 
0.66: 
0.66: 
0.47: 
0.50: 
0.14: 
0.82: 
0.78: 
0.61: 
0.34: 



0.81 
0.15 
0.15 
0.31 
0.56 
1.41 
0.16 
0.12 
0.22 
1.09 



Reference 



References. (l) |Sousaet aLlpOTT) ; (2) |Santos et"aT] ( |2004) ; (3) |Segransan et aL|p0TT] >; (4) |Santos etaT7| ( [2005] >. 



0.4 



•3 0.2 



-0.2 





- 



4400 



4600 



4800 



5000 



T e „ (K) 



Fig. 9. [Cr i/Cr n] and [Ti i/Ti n] as a function of effective tem- 
perature. The black squares and red asterisks correspond to the 
abundance ratios derived using old and new stellar parameters, 
respectively. The solid and dashed lines depict the linear fits of 
the data. 



-0.16 ± 0.08, whereas the slope with the parameters of SO08 
is -0.49 ± 0.08. The new slope of [Til/Tin] is also improved 
a lot and is -0.18 + 0.06 dex per 1000 K. For comparison the 
slope of SO08 is -0.39 ± 0.08 dex. Although the trend with T eS 
is weak, there is a shift of about 0.2 dex for the [Ti i/Ti n] ratio. 
This shift is difficult to connect to the still possible uncertainties 
in the stellar parameters and its exact nature still remains to be 
clarified. Probably, one (or more) of the above mentioned effects 
can be responsible for that. Unfortunately, in the literature there 
is no available NLTE calculations for the Ti i and Ti n lines used 
in our study, and it is difficult to estimate the NLTE effect for the 
[Tii/Tin] ratio. 

Summarizing, this independent test shows that the new stel- 
lar parameters derived from the iron lines more carefully cho- 
sen for cooler stars make the observed [Cr i/Cr n] and [Ti i/Ti n] 
trends with T e ff much weaker. 



8. Conclusions 

In this work, we present a new iron line list in order to correct 
for discrepancies in temperatures for a sample of 451 stars that is 
part of the HARPS GTO program using the cool star HD 21749 
as a reference. The quality of the line list plays a key role for 
temperature determination especially for the K type stars. The 
new line list is compiled in order to eliminate blended lines in 
the spectra of these stars that suffer more from such effects. We 
also apply a limitation to very strong and very weak lines that 
usually cause errors in the EW measurements (see Sect. [XT). 

We derived the stellar parameters for the 451 stars of the 
sample in a homogeneous way with the new line list and com- 
pare our results with the work of SO08, where the authors fol- 
lowed the same analysis but with an expanded line list. 

We find very good agreement for the high and intermediate 
ranges in temperature with SO08. The differences appear in the 
lower temperatures below ~5000 K. In addition, surface gravity 
and metallicitity remain unaffected with the new line list. These 
results are very important since accurate parameters of planet 
host stars are essential to characterize the planets. In particu- 
lar, the metallicities of this sample are used for the investigation 
between the stars and their planet hosts. The fact that the new 
metallicities are not different from the old ones supports the al- 
ready established correlation between the planet host stars and 
planet frequency. 

The agreement with the IRFM that is considered to be a less 
model dependent technique, suggests that our parameters are 
more precise at least for the cooler stars compared to the previ- 
ous study of SO08 where the existing trend in low temperatures 
between the spectroscopic and the photometric temperatures dis- 
appears with the new line list. In addition, we compare our new 
temperatures with the ones derived with interferometry for 9 star 
with available angular diameter measurements. Even though the 
sample is very small, a comparison between these results shows 
very good agreement. Attention should be made to temperatures 
derived with the angular diameters that have larger uncertainties. 

Having checked the reliability of the new line list, we apply 
it to determine fundamental parameters for some of the planet 
hosts with low temperatures, which were previously analyzed 
with different line lists with available high S/N spectra. Finally, 
we use the new parameters for stars with T e ff <5000 K to sig- 
nificant reduce the trends of [Cr i/Cr n] and [Ti i/Ti n] abundance 
ratios with temperature. 
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Fig. A.l. Reduced EW versus excitation potential for the line list 
of this work (circles) and for SO08 (crosses). The slope of the 
linear fits is also depicted. 
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Appendix A: Reduced EW versus Excitation 
Potential 

We show the correlation between the reduced EW and excitation 
potential (Fig. 



A.l i. For the line list of this work there is much 



smaller correlation than in the work of SO08 suggesting small 
degeneracies between f, and T e g. 



Appendix B: The microturbulence relationship 

Microturbulence is taken into consideration for abundance anal- 
yses to reconcile differences between the observed and predicted 
from models equivalent widths of strong lines. Previous studies 
of FGK dwarfs have shown that f, depends on T e g and logg 
(e.g.|Nissen|198 1 ; Reddy et al. 2003 ; Allen de Prieto et al.|2004 



Adibekyan et al. 2012a; Ram irez et al.|2013 i. Using a linear re- 
gression analysis to the new parameters of the sample, we derive 
the following expression: 



ft = 6.932 x 10" 4 r eff - 0.348 logg - 1.437. 



(B.l) 
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Fig. B.l. Correlation of microturbulence with temperature and 
surface gravity as described by equation B.l We set log g to the 
highest, average and lowest values of the sample. 



Here, £, is in km/s, T e ff and log g are in their traditional units. 
The parameters of the stars in the sample range: 

4400 < r eff < 6400 K, 

3.6 < logg <4.8 dex, 

-0.8 < [Fe/H] < 0.4 dex. 

The new derived parameters indicate a linear dependence on 
temperature for a set value of surface gravity. In Fig. B.l we see 



the dependence of microturbulence on temperature for a set of 
logg values. Microturbulence clearly increases with temperature 
and decreases with surface gravity. 
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